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Abstract. Due to the isotropy of d-dimensional hyperspherical space, one expects there to 
exist a spherically symmetric fundamental solution for its corresponding Laplace-Beltrami 
operator. The i?-radius hypersphere with i? > 0, represents a Riemannian manifold 
with positive-constant sectional curvature. We obtain a spherically symmetric fundamental 
solution of Laplace's equation on this manifold in terms of its geodesic radius. We give 
several matching expressions for this fundamental solution including a definite integral over 
reciprocal powers of the trigonometric sine, finite summation expressions over trigonometric 
functions. Gauss hypergeometric functions, and in terms of the associated Legendre function 
of the second kind on the cut (Ferrers function of the second kind) with degree and order 
given by d/2 — 1 and 1 — d/2 respectively, with real argument between plus and minus one. 

Key words: hyperspherical geometry; fundamental solution; Laplace's equation; separation 
of variables; Ferrers functions 
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1 Introduction 

We compute closed-form expressions of a spherically symmetric Green's function (fundamental 
solution of Laplace's equation) for a d-dimensional Riemannian manifold of positive-constant 
sectional curvature, namely the i2-radius hypersphere with i? > 0. This problem is intimately 
related to the solution of the Poisson equation on this manifold and the study of spherical har- 
monics which play an important role in exploring collective motion of many-particle systems 
in quantum mechanics, particularly nuclei, atoms and molecules. In these systems, the hy- 
perradius is constructed from appropriately mass-weighted quadratic forms from the Cartesian 
coordinates of the particles. One then seeks either to identify discrete forms of motion which 
occur primarily in the hyperradial coordinate, or alternatively to construct complete basis sets 
on the hypersphere. This representation was introduced in quantum mechanics by Zernike &: 
Brinkman j37j . and later invoked to greater effect in nuclear and atomic physics, respectively, 
by Delves [5] and Smith |31] . The relevance of this representation to two-electron excited states 
of the helium atom was noted by Cooper, Fano & Prats [3]; Fock [lOj had previously shown that 
the hyperspherical representation was particularly efficient in representing the helium wave- 
function in the vicinity of small hyperradii. There has been a rich literature of applications 
ever since. Examples include Zhukov [38] (nuclear structure), Fano [9] and Lin j24j (atomic 
structure), and Pack &: Parker |28] (molecular collisions). A recent monograph by Berakdar [2] 
discusses hyperspherical harmonic methods in the general context of highly-excited electronic 
systems. Useful background material relevant for the mathematical aspects of this paper can be 
found in [22l [331 [M] • Some historical references on this topic include [171 (SS [291 [SOI [36] • 
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This paper is organized as follows. In Section [2] we describe hyperspherical geometry and its 
corresponding metric, global geodesic distance function, Laplacian and hyperspherical coordinate 
systems which parametrize points on this manifold. In Section [3] for hyperspherical geometry, 
we show how to compute 'radial' harmonics in a geodesic polar coordinate system and derive 
several alternative expressions for a 'radial' fundamental solution of the Laplace's equation on 
the i?-radius hypersphere. 

Throughout this paper we rely on the following definitions. For ai, 02, . . . £ C, i,j G Z and 
j j 

j < i then ^ a„ = and ]^ a„ = 1. The set of natural numbers is given by N := {1, 2, 3, . . .}, 
the set No := {0, 1, 2, . . .} = N U {0}, and the set Z := {0, ±1, ±2, . . .}. 

2 Hyperspherical geometry 

The Euclidean inner product for R'^"''^ is given by (x, y) = xoyo + xiyi + • • • + XdUd- The variety 

(x, x) = XQ + xf-l h = i?^, for X G R'^^^ and R> 0, defines the i?-radius hypersphere S^. 

We denote the unit radius hypersphere by S'^ := Sf. Hyperspherical space in d-dimensions, de- 
noted by S^, is a maximally symmetric, simply connected, d-dimensional Riemannian manifold 
with positive-constant sectional curvature (given by l/R'^, see for instance [22l p. 148]), whereas 
Euclidean space R'^ equipped with the Pythagorean norm, is a Riemannian manifold with zero 
sectional curvature. 

Points on the d-dimensional hypersphere can be parametrized using subgroup-type coor- 
dinate systems, i.e., those which correspond to a maximal subgroup chain 0{d) D • • • (see for 
instance [lEl [20]). The isometry group of the space is the orthogonal group 0(d). Hyper- 
spherical space S^, can be identified with the quotient space 0{d)/0{d — 1). The isometry 
group 0{d) acts transitively on S^. There exist separable coordinate systems on the hyper- 
sphere, analogous to parabolic coordinates in Euclidean space, which can not be constructed 
using maximal subgroup chains. Poly spherical coordinates, are coordinates which correspond to 
the maximal subgroup chain given by 0(d) D • • • . What we will refer to as standard hyper- 
spherical coordinates, correspond to the subgroup chain given by 0{d) D 0{d— 1) D • • • D 0(2). 
(For a thorough discussion of polyspherical coordinates see Section IX. 5 in [35] •) Polyspherical 
coordinates on all share the property that they are described by (d+ l)-variables: R £ [0, 00) 
plus d-angles each being given by the values [0,27r), [0,7r], [— 7r/2,7r/2] or [0, 7r/2] (see [TSl [T9]). 

In our context, a useful subset of polyspherical coordinate are geodesic polar coordinates 
(^,x) (see for instance [27]). These coordinates, which parametrize points on S^, have origin 
at O = (i?, 0, . . . , 0) G R'^+^ and are given by a 'radial' parameter 9 G [0, vr] which parametrizes 
points along a geodesic curve emanating from O in a direction x G S*^"^. Geodesic polar 
coordinate systems partition into a family of {d — l)-dimensional hyperspheres, each with 
a 'radius' := 6d & (0, vr), on which all possible hyperspherical coordinate systems for S'^"^ may 
be used (see for instance [35]). One then must also consider the limiting case for = 0, vr to 
fill out all of S^. Standard hyperspherical coordinates (see \21\ 125]) are an example of geodesic 
polar coordinates, and are given by 

xq = Rcos 9, 

xi = Rsm9 cos6d-i, 

X2 = R sin 9 sin 9d-i cos 9d-2, 

(1) 

Xd-2 = R sin 6 sin 9d~i ■ ■ ■ cos 92 , 
Xd~i = R sin 6 sin 9d~i ■ ■ ■ sin 92 cos 0, 
Xd = R sin sin 9d-i ■ ■ ■ sin 92 sin (f). 



9i G [0, vr] for i G {2, . . . , d], 9 := 9d, and (j) G [0, 2-k). 
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In order to study a fundamental solution of Laplace's equation on the hypersphere, we need 
to describe how one computes the geodesic distance in this space. Geodesic distances on are 
simply given by arc lengths, angles between two arbitrary vectors, from the origin in the ambient 
Euclidean space (see for instance \22\ p. 82]). Any parametrization of the hypersphere S^, must 
have (x, x) = Xq + ■ ■ ■ + x"^ = R"^, with R > 0. The distance between two points x, x' S on 
the hypersphere is given by 

,M = Ry = Rcos- {j;^^^) = • (2) 

This is evident from the fact that the geodesies on are great circles, i.e., intersections of 
with planes through the origin of the ambient Euclidean space, with constant speed parametriza- 
tions. 

In any geodesic polar coordinate system, the geodesic distance between two points on the 
submanifold is given by 



(i(x, x') = i? cos ^ ^-^(x, x')^ = Rcos ^ (cos 6* cos 0' + sin 6* sin 0' cos 7) , 



(3) 



where 7 is the unique separation angle given in each polyspherical coordinate system used to 
parametrize points on S*^"^. For instance, the separation angle 7 in standard hyperspherical 
coordinates is given through 

d-2 d-2 i~l 

COS 7 = cos((/) — (j)') sin 0jsin 9/ + cos 6'jcos 9/ sin 9jsin 9/. (4) 

i=l i=l j = l 

Corresponding separation angle formulae for any hyperspherical coordinate system used to 
parametrize points on S'^"^ can be computed using ([2]) and the associated formulae for the 
appropriate inner-products. 

One can also compute the Riemannian (volume) measure dvolg (see for instance Section 3.4 
in [E]), invariant under the isometry group SO{d), of the Riemannian manifold S^. For instance, 
in standard hyperspherical coordinates ^ on the volume measure is given by 

dvolg = R"^ sin°'-^ 6 d9 duj := R'^ sin'^-^ 9 d9 sin'^^^ 9d^i • • • sin ^2 d9i--- d9d^i. (5) 

The distance r G [0, 00) along a geodesic, measured from the origin, is given by r = 9R. To show 
that the above volume measure ([5]) reduces to the Euclidean volume measure at small distances 
(see for instance [21]), we examine the limit of zero curvature. In order to do this, we take the 
limit — 7- O"*" and 72 — )• 00 of the volume measure ([5]) which produces 

dvoL ~ R'^-^ sin'^-^ ( —) drdu ~ r^^^drdoj, 



which is the Euclidean measure in R"^, expressed in standard Euclidean hyperspherical coordi- 
nates. This measure is invariant under the Euclidean motion group E{d). 

It will be useful below to express the Dirac delta function on S^. The Dirac delta function on 
the Riemannian manifold with metric g is defined for an open set U C with x, x' € 
such that 



6gMd.Olg=\\ (6) 



For instance, using ([5]) and ([6]), in standard hyperspherical coordinates on ([T]), we see that 
the Dirac delta function is given by 



5{9-9') 5{9i-9[)---6{9d^i-e'^^^, 
R<^ sin'^-i 9' sin 6*^ • • • sin'^^^ ^/^_^ 



5g{x,x') 
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2.1 Laplace's equation on the hypersphere 

Parametrizations of a submanifold embedded in Euclidean space can be given in terms of coor- 
dinate systems whose coordinates are curvilinear. These are coordinates based on some trans- 
formation that converts the standard Cartesian coordinates in the ambient space to a coordinate 
system with the same number of coordinates as the dimension of the submanifold in which the 
coordinate lines are curved. 

The Laplace-Beltrami operator (Laplacian) in curvilinear coordinates = {C^ , ■ ■ ■ , C^) on 
a Riemannian manifold is given by 

1 d r ^ d 



where Ig] = | det {gij)\, the metric is given by 

d 



ds^ = Yl 9^Jdede, (8) 



and 



1=1 

where G {0, 1} is the Kronecker delta 
1 I 1 if * = 7, 

si--=\. (9) 

[0 ifz/j, 

for i,j G Z. The relationship between the metric tensor dj = diag(l,... ,1) in the ambient 
space and gij of ([7|) and dH) is given by 



The Riemannian metric in a geodesic polar coordinate system on the submanifold is 
given by 

ds^ = {d9^ + sin^ 9 d-f^) , (10) 

where an appropriate expression for 7 in a curvilinear coordinate system is given. If one com- 
bines ([I]), (01), dZI) and ([To]) , then in a geodesic polar coordinate system, Laplace's equation on 
is given by 



J +(d-l)cot0-^ + -^As.-i/ 



sm 



0, (11) 



where Agd-i is the corresponding Laplace-Beltrami operator on S*^ ^. 
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3 A Green's function on the hypersphere 
3.1 Harmonics in geodesic polar coordinates 

The harmonics in a geodesic polar coordinate system are given in terms of a 'radial' solution 
('radial' harmonics) multiplied by the angular solution (angular harmonics). 

Using polyspherical coordinates on S"^"^, one can compute the normalized hyperspherical 
harmonics in this space by solving the Laplace equation using separation of variables. This 
results in a general procedure which, for instance, is given explicitly in I19j . These angu- 
lar harmonics are given as general expressions involving trigonometric functions, Gegenbauer 
polynomials and Jacobi polynomials. The angular harmonics are eigenfunctions of the Laplace- 
Beltrami operator on S*^"^ which satisfy the following eigenvalue problem (see for instance (12.4) 
and Corollary 2 to Theorem 10.5 in [32]) 

As.-iy;^(X) = -/(/ + d - 2)y,^(S), (12) 

where x S S*^"^, yj^(x) are normalized angular hyperspherical harmonics, / G Nq is the angular 
momentum quantum number, and K stands for the set of {d — 2)-quantum numbers identifying 
degenerate harmonics for each I and d. The degeneracy 

(see (9.2.11) in [35]), tells you how many linearly independent solutions exist for a particular I 
value and dimension d. The angular hyperspherical harmonics are normalized such that 



/ y,^(X)l^F(X)da; = 5f5 



K ) 



where du) is the Riemannian (volume) measure on S"^"^, which is invariant under the isometry 
group SO{d) (cf. dSj)), and for x + iy = z G C, z = x — iy, represents complex conjugation. 
The angular solutions (hyperspherical harmonics) are well-known (see Chapter IX in j35j and 
Chapter 11 [8]). The generalized Kronecker delta symbol 5^ (cf. ([9|)) is defined such that it 
equals 1 if all of the {d — 2)-quantum numbers identifying degenerate harmonics for each / and d 
coincide, and equals zero otherwise. 

We now focus on 'radial' solutions of Laplace's equation on S^, which satisfy the following 
ordinary differential equation (cf. ([TT]) and p^ ) 



d^u , , , du 1(1 + d -2) , , 

+ (d- l)cot6'— - ^ ^ti = 0. (13) 



de^ ^ ' dO 



sm 



Four solutions of this ordinary differential equation u^!^,U2l. '■ (—1, 1) — )• C are given by 
and 

^2i(cos0) := (^.^gj,/,_, Qgyr'^'^(^°^^)^ (14) 

where PlJ,Q(^ : (—1,1) — t- C are Ferrers functions of the first and second kind (associated 
Legendre functions of the first and second kind on the cut). The Ferrers functions of the first 
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and second kind (see Chapter 14 in [26]) can be defined respectively in terms of a sum over two 
Gauss hyper geometric functions, for all i', fj, £ C such that u + ji ^ — N, 



2^+1 

P^(x) := — — sm 



2(^ + /^) 



2^ 

H ;= cos 



vr 



r 



2 



:(1 - x')->^l\F, 
(1 - 



1 — — ^ — + 2 3 2 
2 ' 2 '2'^ 



-U — U — fi + 1 1 2 

~"2 ' 2 ' 2'^ 



(cf. (14.3.11) in 



and 



"TT 



Qi;(x) := V^2'^cos + 



{1 - x^yf^/^Fi 



1 — 1/ — u — 11 + 2 3 



2'^ 



^/^2■ 



sm 



vr 



2 



(1 - x2)-^/22Fi 



-U — fl U — fl + 1 1 



' 2'' 



(15) 



(cf. (14.3.12) in [26]). The Gauss hypergeometric function 2-F1 : C x C x (C \ -Nq) x {z g C : 
|z| < 1} — )• C, can be defined in terms of the infinite series 

2F1 [a, b; c; z) := — -z 



n=0 



(c)„n! 



(see (15.2.1) in [26]), and elsewhere in z by analytic continuation. On the unit circle \z\ = 1, the 
Gauss hypergeometric series converges absolutely if Re (c — a — 6) € (0, 00), converges conditio- 
nally if z / 1 and Re (c — a — 6) G (—1,0], and diverges if Re (c — o — 6) G (—00, —1]. For z G C 
and n G Nq, the Pochhammer symbol (z)„ (also referred to as the rising factorial) is defined as 
(cf. (5.2.4) in 



On :='[liz + i-l). 



i=l 



The Pochhammer symbol (rising factorial) is expressible in terms of gamma functions as (5.2.5) 
in [26] 



T{z + n) 

TizT'' 



for all z G C \ — Nq. The gamma function F : C \ — Nq — )• C (see Chapter 5 in |26j ) is an 
important combinatoric function and is ubiquitous in special function theory. It is naturally 
defined over the right-half complex plane through Euler's integral (see (5.2.1) in [26]) 



r{z) := / 
^0 



Rez > 0. The Euler reflection formula allows one to obtain values of the gamma function in 
the left-half complex plane (cf. (5.5.3) in [26] ) , namely 



F(z)F(l-z) 



TT 



smvTZ 
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< Rez < 1, for Rez = 0, z / 0, and then for z shifted by integers using the following 
recurrence relation (see (5.5.1) in |26| ) 

r(z + i) = zr(z). 

An important formula which the gamma function satisfies is the duplication formula (i.e., (5.5.5) 
in [26]) 

222-1 / 1\ 

r(2z) = ^r(z)r(z + -j, (i6) 



provided 2z — Nq. 

Due to the fact that the space is homogeneous with respect to its isometry group, the 
orthogonal group 0{d), and therefore an isotropic manifold, we expect that there exist a fun- 
damental solution on this space with spherically symmetric dependence. We specifically expect 
these solutions to be given in terms of associated Legendre functions of the second kind on the 
cut with argument given by cos 9. This associated Legendre function naturally fits our require- 
ments because it is singular at 9 = 0, whereas the associated Legendre functions of the first 
kind, with the same argument, is regular at ^ = 0. We require there to exist a singularity at 
the origin of a fundamental solution of Laplace's equation on S^, since it is a manifold and 
must behave locally like a Euclidean fundamental solution of Laplace's equation which also has 
a singularity at the origin. 

3.2 Fundamental solution of the Laplace's equation on the hypersphere 

In computing a fundamental solution of the Laplacian on S^, we know that 

-A5^(x,x') = 53(x,x'), (17) 

where g is the Riemannian metric on (e.g., (jlOp ) and 6g is the Dirac delta function on the 
manifold (e.g., dS])). In general since we can add any harmonic function to a fundamental 
solution of Laplace's equation and still have a fundamental solution, we will use this freedom 
to make our fundamental solution as simple as possible. It is reasonable to expect that there 
exists a particular spherically symmetric fundamental solution Sf^ on the hypersphere with pure 
'radial', 9 := (i(x, x') (e.g., and constant angular dependence due to the influence of the 
point-like nature of the Dirac delta function in (I17p . For a spherically symmetric solution to 
the Laplace equation, the corresponding Agd-i term in (jlip vanishes since only the I = term 
survives in (|13p . In other words we expect there to exist a fundamental solution of Laplace's 
equation on such that 5^(x, x') = f{9) (cf. where Ris a parameter of this fundamental 
solution. 

We have proven that on the i2-radius hypersphere S^, a Green's function for the Laplace 
operator (fundamental solution of Laplace's equation) can be given as follows. 

Theorem 1. Let d € {2, 3, . . .}. Define Id : (0, it) ^ H as 

dx 



■ fi—i ' 
sm X 



x,x' G S^, and Sf^ : {S% x S^) \ {(x,x) : x G S^} R defined such that 
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where 9 := cos~^ ([x, x']) is the geodesic distance between x and x' on the unit radius hyper- 
sphere S'^, with X = x/i?, x' = x'/i2, then 5^ is a fundamental solution for —A where A is the 
Laplace-Beltrami operator on S^. Moreover, 



{d-2)\l 

d- 



9 /^'(2A:-2)!! 1 ' 

logcot- + cos^g p-I)TT-^ 



if d even, 



(d-l)/2 



cot 



2fc-l , 



or 

(^-3)!! 

{d-iy.i 



(2A;- 1)(A:- l)!(((i-2A:- l)/2) 

(d-l)/2 



COS 9 



(2A;-3)!! 



1 



fc=i 



(2fc-2)!!sin2'=-i0' 



if d odd, 



„ ^ , 1 d 3 2 
cos 6' 2-^1 ( 2' 2' 2'^°*^ ' 



cos I 



sm 



^ , 3 — d 3 2 
2^1 1, — — ; -;cos I 



, r (d/2) 2''/2-i (sin0)'^/2-i ^'^/2-i^^°^ '^^ 

In the rest of tliis section, we develop the material in order to prove this theorem. 
Since a spherically symmetric choice for a fmidamental solution satisfies Laplace's equation 
everywhere except at the origin, we may first set g = f in (|lip and solve the first-order equation 

g' + {d- 1) cos9 g = 0, 

which is integrable and clearly has the general solution 



9(0) 



d9 



Co (sin ( 



(18) 



where cq € R is a constant. Now we integrate (jlSp to obtain a fundamental solution for the 
Laplacian on 



5^(x,x') = coXd(0) + ci, 
where Xd : (0, vr) — )■ R is defined as 

™/2 dx 



MO) 



sm X 



(19) 



(20) 



and co,ci S R are constants which depend on d and R. Notice that we can add any harmonic 
function to (jl9p and still have a fundamental solution of the Laplacian since a fundamental 
solution of the Laplacian must satisfy 



(-A<^)(x')5^(x,xOdvol' =(/p(x). 



for all if € S{Sf,), where 5 is the space of test functions, and dvol^ is the Riemannian (volume) 
measure on in the primed coordinates. Notice that our fundamental solution of Laplace's 
equation on the hypersphere ( (fT9|) . (pO|) ) has the property that it tends towards +oo as 6* ^ 0+ 
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and tends towards — oo as — t- vr". Therefore our fundamental solution attains all real values. 
As an aside, by the definition therein (see |14[ I15j). is a parabolic manifold. Since the 
hypersphere is bi-hemispheric, we expect that a fundamental solution of Laplace's equation 
on the hypersphere should vanish at 6 = 7r/2. It is therefore convenient to set ci = leaving us 
with 

5^(x,x') = coXdW. (21) 

In Euclidean space R^, a Green's function for Laplace's equation (fundamental solution for 
the Laplacian) is well-known and is given by the following expression (see p. 94], [12\ p. 17], 
[31 p. 211], [61 p. 6]). Let d G N. Define 



a''(x,x')= < 



r(d/2) „„ ^/||2-d ifd=iord>3, 



2^^/2((i-2)" " - ' ^22) 

— logllx-x'ir^ ifd = 2, 



then Q'^ is a fundamental solution for —A in Euclidean space R*^, where A is the Laplace 
operator in R'^. Note that most authors only present the above theorem for the case d> 2 but 
it is easily-verified to also be valid for the case d = 1 as well. 

The hypersphere S^, being a manifold, must behave locally like Euclidean space R"'. There- 
fore for small 9 we have ~ 1 + and ~ 1 — and in that limiting regime 



-log 6* ifd = 2, 
^ if d > 3, 



which has exactly the same singularity as a Euclidean fundamental solution. Therefore the 
proportionality constant cq is obtained by matching locally to a Euclidean fundamental solution 

5^ = cold - (23) 

in a small neighborhood of the singularity at x = x', as the curvature vanishes, i.e., R — )• oo. 

We have shown how to compute a fundamental solution of the Laplace-Beltrami operator 
on the hypersphere in terms of an improper integral (I20p . We would now like to express this 
integral in terms of well-known special functions. A fundamental solution can be computed 
using elementary methods through its definition ()20p . In d = 2 we have 

r/2 dx 1 cose + 1 e 

M^) = / - — = t; log n — T = o' 

Jq sm X 2 cos e — 1 2 

and in d = 3 we have 

i3{e)= / -V = cote. 

Je sin^ X 
In d £ {4, 5, 6, 7} we have 

^ /^\ 1 1 ^ cos^ 

X4(6') = - log cot - + 



2 ^ 2 2sm^e' 

1^(9) = cot6l-F -cot^6l, 
3 

_ 3- e 3cos0 cose 

laie) = - log cot - + . a + . and 
o 2 8 sm e 4 sm e 

2 1 

1^(0) = cote + - cot^ e + - cot^ e. 

3 5 
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Now we prove several equivalent finite summation expressions for ld{0). We wish to compute 
the antiderivative Jm '■ (0, vr) — t- R, which is defined as 

dx 



sin'"^' 



where m G N. This antiderivative satisfies the following recurrence relation 

cos X 



■^m (^) 



+ 1 TT^m-2{x), 



(m — l)sin'^ ^x ' (m — 1) 
which follows from the identity 



(24) 



+ 



cos x 



sin"^x sin"* ' sin™x 



cos X, 



and integration by parts. The antiderivative 3m{x) naturally breaks into two separate classes, 
namely 



dx 



and 



dx 



{2n - 1) 
(2n)!! 

' (2n-2)!! 
'(2n - 1) 



X ^(2k-2)\l 1 

log cot — h cos X > — — ^7 

^ 2 i2k - ly.l sin^'' X 



+ C, 



(25) 



sm^" X 



COS X 



E 



(2fc-3)!! 



n 



(2A;-2)!!sin2'=-ix 
cot^'^"^ X 



or 



(26) 



k=l 



{2k - l){k - iy.{n - k)\ 



+ c, 



where C is a constant. The double factorial (•)!! : { — 1, 0, 1, ...}—> N is defined by 

n • (n — 2) • • • 2 if n even > 2, 
n • (n — 2) • • • 1 if n odd > 1, 
1 ifnG{-l,0}. 



n! 



II ■- 



Note that (2n)!! = 2"'n! for n € Nq. The finite summation formulae for 'Jmix) all follow trivially 
by induction using (j24p and the binomial expansion (cf. (1.2.2) in j26j ) 



2 



(1 + COS X 



n 



k=0 



cot^*^ X 
kl{n-k)\' 



The formulae ([25]) and ([26|) are essentially equivalent to (2.515.1-2) in [13], except (2.515.2) is 
in error with the factor 28^^ being replaced with 2^. This is also verified in the original citing 
reference [3l]. By applying the limits of integration from the definition of X^(^) in (j20p to (j25p 
and ()26p we obtain the following finite summation expression 



(d-2)!! 



d/2~l 



{2k-2)\\ 1 



log cot — \-cos9 , , „, 
^2 ^ (2fc- l)!!sin2*^0 



if d even. 



d-3 



or 

(d-3)!! 



(d-l)/2 



cot 



2/c-l 



^ (2A;- 1)(A;- l)!((d-2A;- l)/2)!^ 



(27) 



> if d odd. 



(d-l)/2 



(2A;-3)!! 



1 



(^-2)!!^^°'^ (2fc-2)!!sin2^-i0' 
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Moreover, the antiderivative (indefinite integral) can be given in terms of the Gauss hyper- 
geometric function as 

^ = - cos 2^1 a ^; I, cos' e) + C, (28) 



sm 



2' 2' 2' 



where C G R. This is verified as follows. By using 

2Fi(a, b; c; z) = — 2-Pi(a + 1, 6 + 1; c + 1; z) 
az c 

(see (15.5.1) in [2B]), and the chain rule, we can show that 
— cos t) 2-ri cos 



d9 V2 2' 2 

,\ ^ fl d 3 .\ d 0, ^ [3 d + 2 5 n 
= sm6' 2^1 I 2' 2' 2''^°'^ I + - cos 2-F1 I -, ; -; cos 

The second hypergeometric function can be simplified using Gauss' relations for contiguous 
hypergeometric functions, namely 

z 2Fi{a + l,b+l;c+l;z) = -^\2Fi{a,b + 1; c; z) - 2^1(0 + 1, 6; c; z)l 
(see O p. 58]), and 

5 — a ^ , , , a 



2 



Fi(a,6 + l;c; z) = ——2Fi{a,b;c;z) + - 2-^1(0 + l-,b;c;z) 



(see (15.5.12) in [5B]). By applying these formulae, the term with the hypergeometric function 
cancels leaving only a term which is proportional to a binomial through 

2F,{a,b;b;z) = {l-z)-'' 

(see (15.4.6) in [26]), which reduces to l/sin"^"^ 9. By applying the limits of integration from the 
definition of X(i{9) in (|20|) to (|28p we obtain the following Gauss hypergeometric representation 



Xrf(0)=COSe2i^lQ,^;|cOs2ey 



(29) 



Using (j29p . we can write another expression for Id{9). Applying Eulers's transformation 

2F1 (a, 6; c; z) = (1 - z)""""^ 2^^! (c - a, c - 6; c; z) 
(see (2.2.7) in [1]), to ([29]) produces 

cos „ / 3 — d 3 9 



= . rf_o^ 2i^i 1, -; cos 



sm'^ ^ 



2 ' 2' 



Our derivation for a fundamental solution of Laplace's equation on the i2-radius hypersphere 
in terms of Ferrers function of the second kind (associated Legendre function of the second kind 
on the cut) is as follows. If we let v + = in the definition of Ferrers function of the second 
kind Q'i : (-1, 1) ^ C 1^, we derive 



Q-/{x) = ^ 'yt- 2Fi \ -,v + X' 
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for all G C. If we let = d/2 — 1 and substitute x = cos 9, then we have 



IT cosfc'sm'^/^-ig) _ fl d 3 2 



Q,-"^^(cos0) = ^;/2-i r (^) U' 2' 2'"°^ V • 

Using the duplication formula for gamma functions ()16p . then through pOD we have 

We have therefore verified that the harmonics computed in Section [3.H namely (O) 
an alternate form for a fundamental solution of the Laplacian on the hypersphere. 

Note that as a result of our proof, we see that the relevant associated Legendre functions of 
the second kind on the cut for d G {2, 3, 4, 5, 6, 7} are (cf. (j27|) and ([3T]) ) 

Q 

Qo(cos6') = log cot -, 



(sin ^) 1/2 ^1/2 v--'^;- Y 2 

1 / ^^ 1 , ^ COS 

-Qi (cos 6') = -logcot- + 



sine 2 """""^ 2 ■ 2sin2e' 

^ -Q"f/^(cos 61) = ^ W J (cot 6* + ^ cot^ 



(sin 0)3/2 "^3/2 '^^^"-^ 2 V 2 V 3 

1 ^„2/ /IN 1 T ^ COS0 COS0 , 

rrrrQ^ ( COS = - log cot — I o 1 ^ — , and 

(sin 0)2 ^2 V ^8^2 8sin20 12sin^0 



^ 5/2/_ /)N _ 1 /^_. /) , 2 3 1 5 



(cos0) = - W- COt0 + - cot30 + - COt^0 



(sin 0)5/2 ^5/2 sY 2 (^--^--3--^ '^-5 

The constant cq in a fundamental solution for the Laplace operator on the hypersphere (j2ip 
is computed by locally matching up, through (|23p . to the singularity of a fundamental solution 
for the Laplace operator in Euclidean space (122]) . The coefficient cq depends on d and -R. For 
d > 3 we take the asymptotic expansion for cqT(i{0) as — )• 0"*", and match this to a fundamental 
solution for Euclidean space (f22]l . This yields 

r(^/2) ,„^, 

For d = 2 we take the asymptotic expansion for 

£0X2(6') = -cologtan- ~ colog||x-x ||~\ 

as —7- 0"*", and match this to ^^(x, x') = (2ir)~^ log ||x — x'||~^, therefore cq = {2tt)~^ . This 
exactly matches (j32p for d = 2. The R dependence of cq originates from ()20p . where x and 
represents geodesic distances (cf. ([3])). The distance r S [0,oo) along a geodesic, as measured 
from the origin of S^, is given by r = OR. To show that a fundamental solution (j2ip reduces 
to the Euclidean fundamental solution at small distances (see for instance [21]), we examine 
the limit of zero curvature. In order to do this, we take the limit 0—7-0"'' and i? — t- 00 of (j20p 
with the substitution x = r/R which produces a factor of R'^~^. So a fundamental solution of 
Laplace's equation on the Riemannian manifold is given by 

^<i(^ - r(rf/2) ^ 

The proof of Theorem [1] is complete. 
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Apart from the well-known historical results in two and three dimensions, the closed form 
expressions for a fundamental solution of Laplace's equation on the i?-radius hypersphere given 
by Theorem [1] in Section 13.21 appear to be new. Furthermore, the Ferrers function (associated 
Legendre) representations in Section 13.11 for the radial harmonics on the i?-radius hypersphere 
do not appear to be have previously appeared in the literature. 
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